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ABSTRACT 

The spectra of high-energy protons and nuclei accelerated by supernova remnant shocks are 
calculated taking into account magnetic held amplihcation and Alfvenic drift both upstream and 
downstream of the shock for different types of supernova remnants during their evolution. The 
maximum energy of accelerated particles may reach 5 • 10 18 eV for Fe ions in Type lib SNRs. 
The calculated energy spectrum of cosmic rays after propagation through the Galaxy is in good 
agreement with the spectrum measured at the Earth. 

Subject headings: acceleration of particles — shock waves — supernova remnants 



1. Introduction 

Supernova remnants (SNRs) are recognized as 
the principal sources of Galactic cosmic rays, and 
the diffusive shock acceleration is accepted as the 
mechanism of cosmic ray acceleration by a su- 
pernova blast wave moving through the turbulent 
interstellar medium. Accelerated by supernova 
shocks, the energetic particles diffuse in the inter- 
stellar magnetic fields and fill the entire Galaxy. 
Clear evidence for particle acceleration in SNRs is 
given by observations of non-thermal radio, X-ray, 
and gamma-ray radiation. It is experimentally 
established from H.E.S.S. (High Energy Stereo- 
scopic System) data that there are cosmic-ray 
particles with energies exceeding 10 14 eV in the 
shell of the supernova remnant RX JI7I3. 7-3946 
(jAharonian et al.|[2007h ■ 

For a time, the theoretical estimates of maxi- 
mum pro ton energy were at a level o f -E m ax~10 13 ~ 
10 14 eV (lLagage k Cesarsky|[l983l) . The work of 
Bell ( 2004h demonstrated that a very large ran- 



dom magnetic field SB ^> B lsm (where B{ s - m Ki'b^,G 
is the average interstellar field upstream of the 
shock) can be generated by a cosmic-ray stream- 
ing instability in the precursor of strong shocks. 
It results then in efficient confinement of energetic 
particles in the shock vicinity, and it may raise the 
maximum particle energy by about two orders of 
magnitude (the exact value depends on the super- 
nova parameters, see below). Remarkably, data 
on synchrotron X-ray emission from a number of 
young SNRs proved the pr esence of strong mag- 
netic fields 150 to 500 /xG jVolk et al.l 120051 ) that 
can be naturally explained by the effect of cosmic- 
ray streaming instability. An excellent concise re- 
view on the origin of observ e d cos mic-ray spec- 
trum was presented bv lHillasI (|2006l ). 

In the present work we calculate the steady 
state spectrum of cosmic rays in the Galaxy us- 
ing recent results on magnetic field amplification 
in SNRs and including different types of SNRs in 
the consideration. A numerical code is employed 
in the simulations of particle acceleration and the 
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SNR shock evolution. 

2. Modelling of cosmic ray acceleration in 
supernova remnants 

Because of high efficiency of shock acceleration, 
the spectrum of cosmic rays should be selfcon- 
sistently determined with the account of shock 
modification caused by the pressure of acceler- 
ated particles. We are studing cosmic ray ac- 
celeration and the evolution of a supernova blast 
wave with the use of our numerical c ode de- 
scribed in (| Zirakashvili fc Ptuskinl l2009alibh . The 
hydrodynamic equations are solved together with 
the diffusion-convection transport equation for the 
cosmic ray distribution function f(t,r,p), which 
depends on time t, radial distance from the point 
of supernova explosion r (here the spherical sym- 
metry is assumed), and the particle momentum p. 
The full system of equations is the following: 
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Here p is the gas density, u is the gas velocity, P g 
is the gas pressure, P c = An J p 2 dpvpf/3 is the 
cosmic ray pressure, w(r, t) is the advective veloc- 
ity of cosmic rays, 7 ff is the adiabatic index of the 
gas, and D(r,t,p) is the cosmic ray diffusion coef- 
ficient. It is assumed that the diffusive streaming 
of cosmic rays results in the generation of mag- 
netohydrodynamic waves and provides the Bohm 
diffusion coefficient Db = vpc/ (3ZeB) for accel- 
erating particles of charge Ze and velocity v (c is 
the speed of light). 

The essential feature of our calculations is the 
inclusion of the Alfvenic drift effect for particle 



transport. The Alfven velocity Va = B/^/Anp is 
not negligible in comparison to the gas velocity u 
downstream of the shock if the magnetic field is 
significantly amplified as indicated by the obser- 
vations of the synchrotron X-ray radiation. The 
effective scattering of energetic particles near the 
shock is provided by Alfven waves generated by a 
resonant cosmic ray streaming instability. Unsta- 
ble Alfven waves propagate in the direction oppo- 
site to the cosmic ray gradient. We also assume 
that the growth of waves is balanced by a non- 
linear damping that is strong enough for high wave 
amplitudes. The energy which cosmic rays lose 
due to the wave generation eventually goes into the 
gas heating (see the last term in Eq. (3)). The cos- 
mic ray gradient is negative upstream of the shock 
and the waves generated propagate in the positive 
direction. The waves transmitted downstream of 
the shock are damped and regenerated by the cos- 
mic ray gradient in this region. The cosmic ray 
gradient is positive here because of the adiabatic 
losses of cos mic ray particles (see numeric al cal- 
culations of Zira kashvili fc Ptuskinl (|2009ah . As 
a result the generated Alfven waves propagate in 
the negative direction downstream of the shock. 

The situation is somewhat different for the 
highest energy particles. An electric current of 
these particles amplifies the magnetic field up- 
stream of the sho c k via a non-resonant stream- 
ing instability iBelll (|2004l ). The generated random 
magnetic fields are not Alfven waves but almost 
purely growing magnetic disturbances. This am- 
plified field plays a role of the regular field for low- 
energy particles. They produce Alfven waves that 
propagate along the tangled magnetic field lines. 
The amplified magnetic field is almost isotropic 
even downstream of the shock in spite of compres- 
sion at the shock front. This is because the shock 
front is corrugated by density fluctuations gener- 
ated b y the streaming instability u pstream of the 
shock ( Zirakashvili fc Ptuskinl[2008 ). So the trans- 
port velocity for highest energy particles is close 
to the plasma velocity u in contrast to low-energy 
particles. For the sake of simplicity we neglect this 
below and use the Alfven transport for all energies. 

Taking effects mentioned above into account 
we set the cosmic ray advection velocity equal to 
w = u — Va/V3 downstream of the shock. Because 
of this the accelerating particles "feel" a smaller 
compression ratio and acquire a softer energy spec- 
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trum compared to the usual assumption w = u. 
We set w = u + Va/V3 upstream of the shock 
to account for the Alfvenic drift effect there. The 
damping of Alfven waves created by the cosmic 
ray streaming instability results in a very impor- 
tant effect of gas heating upstream of the shock 
described by the last term in Eq. (3). This effect 
limits the total shock compression ratio. 

The spatial dependence of amplified magnetic 

field is taken in the form B(r) = ^/Airpo m^ Pq 
where po is the gas density of the circumstellar 
medium and Ma is s ome constant. We employ 
results of lVolk et all (j2005h in the analysis of X- 
ray radiation from young SNRs and assume that 
magnetic energy density B 2 /8ir downstream of the 
shock is 3.5 % of the ram pressure pu 2 h that deter- 
mines the constant Ma = 23. It is worth noting 
that this relation is in good agreement with the 
modelling of cosmic ray streaming instab ility in 
young SNRs (|Zirakashvili fc Ptuskinl [20081) . 

The last term in Eq. (4) corresponds to the in- 
jection of thermal protons with momenta p = pi n j 
and mass m at the shock at r = R(t). The di- 
mensionless parameter 77 determines the injection 
efficiency The injection efficiency of thermal ions 
in the process of shock acceleration 7 7 = 0.1u R h/ c 



determined at the time when the Sedov stage be- 
gins. It gives approximately 



is taken in accordance with the paper IZirakashvili 



(|2007h : lere u s h(t) is the time- varying shock ve- 
locity. 

The maximum particle momentum p max reached 
in a process of diffusive shock acceleration can be 
roughly estimated from the condition Dg(p niax ) ~ 
0.1u s i 1 i? s i 1 with Db calculated for the upstream 
magnetic field, which is about 5 times smaller than 
the downstream field. This gives an order of mag- 
nitude estimate p max c/Z ~ 20(u s h/10 3 km s _1 ) 2 i? s h 
TeV, where the shock radius is i? s h pc, and the 
interstellar gas number density n cm~ 3 . 

It c an be shown, see e.g. IPtuskin fc Zirakashvili 
(|2005l ). that the transformation of supernova ex- 
plosion energy to cosmic rays becomes efficient 
from the beginning of the Sedov stage of the shock 
evolution (i.e. when the mass of supernova ejecta 
becomes equal to the mass of swept-up gas) and 
continues later on. As a result, the characteristic 
knee arises in the overall spectrum of particles ac- 
celerated by the evolving supernova remnant. The 
position of knee pknoc can be estimated from the 
above equations for p max where w s h and i? s h are 



Pkn 



3 c/Z ~ 1 • 10 15 £ 5 in 1/6 M e 7 2/3 eV. 



(5) 



Here -E51 is the kinetic energy of the supernova 
explosion in units of 10 51 erg and M e j is the mass 
of supernova ejecta measured in solar masses. 

If the presupernova had a dense star wind with 
velocity u w and the mass loss rate M before the 
explosion, the shock may enter the Sedov stage 
while propagating through the wind material with 
mass density p w = M / {4,-Ku w r 2 ) . Eq. (5) should 
be replaced in this case by the following equation: 



(6) 



fWc/Z ~ 8 • 10 15 S 5 i V /M_ 5 /u U) , 6 MT 1 eV. 

We have performed numerical simulations of 
cosmic ray acceleration for 4 types of supernova 
remnants (they constitute about 90 percents of all 
supernovae). 

1. Type la SNRs with the following parame- 
ters: kinetic energy of explosion E — 10 51 erg, 
number density of the surrounding interstellar gas 
n = 0.1 cm~ 3 , and mass of ejecta M e j = 1.4M©. 
Also important for accurate calculations is the in- 
dex fc, which describes the power law density pro- 



file p s 



of the outer part of the star that 



8M©, and k = 12. 
= 10 51 erg exploding 



freely expands after supernova explosions; k = 7 
for Type la supernova. 

2. Type IIP SNRs with parameters E = 
10 51 erg, n = 0.1 cm~ 3 , M ej = 

3. Type Ib/c SNRs with E 
into the low density bubble with density n = 0.01 
cm~ 3 formed by a progenitor star that starts off 
as an O star, goes through a RSG (Red Super 
Giant) phase , and ends its lif e as a Wolf-Rayet 
star, see e.g. Dwarkad as (2007). The ejecta mass 
is M e j = 2M Q and k = 7. 

4. Type lib SNRs with E = 3 • 10 51 erg, n = 
0.01cm -3 and M e j = 1M Q . Before entering the 
rarefied bubble, the blast wave goes through the 
dense wind emitted by the progenitor star during 
its final RSG stage of evolution. We assume that 
the mass loss rate by the wind is M = 10 -4 MQ/yr 
and the outer wind radius is 5 pc. 

A discussion about properties of SNRs pro- 
d uced by co r e coll apse supernovae can be found 
in lChevalierl (|2005l ). 
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We accept the following relative rates for 
the four types of supernovae described above: 
0.32,0.44,0.22 and 0.02 respectivel y. The first 



three rates are taken from the work of lSmartt et al. 
(2009) based on the statistics of supernovae within 
28 Mpc of the Galaxy (they also ag ree with 
the G alactic supernova rates derived by iLeamanl 
The statistics for rare Type lib events 
is not sufficiently reliable. These supernovae de- 
termine the cosmic ray intensity at energies above 
about 3 x 10 17 eV and our rate 0.02 is chosen to fit 
the data on t hese ultra-high ener gy cosmic rays. 
According to ISmartt et ail (|2009h . the Type lib 
supernovae rate can be as high as 0.04. Notice 
also that the fluctuation effect is very strong for 
these rare events. 

The calculated cosmic ray spectra produced 
over the lifetime of each type of supernovae are 
shown in Figure 1 under the assumption that only 
protons are accelerated. Here Q{p) — 4:np 2 F(p), 
where F(p) is the distribution of all accelerated 
particles injected in the interstellar medium over 
the SNR lifetime. The total number of acceler- 
ated particles is J Q{p)dp. It was assumed that 
the acceleration ceased at t c — 10 5 yr. The shock 
velocity at this moment is close to 200 km/s, and 
the maximum energy of protons confined in the 
supernova remnant is ~ 5 TeV. All particles with 
higher energies accelerated earlier have left the 
remnant. The maximum particle energy at late 
stages of shock evolution can be two-three orders 
of magnitude smaller if one takes into account the 
possible damping of turbulence in the shock pre- 
cursor due to ion-neutral collisions or non-linear 



wave interactions (jPtuskin fc Zirakashvilill2003l ). 

The function F(p) was calculated as the sum of 
two integrals: the integral taken at t c over the vol- 
ume of supernova remnant (47r) JJ^* f(t Cl r,p)r 2 dr, 
and the integral over time of the diffusion flux of 
accelerated particles through the boundary of the 
calculation domain [4ir r 2 J Q C (—Ddf /dr)dt]\ rb . 
The source function Q(p) should be multiplied 
by z/ sn , where v sn is the supernova rate per unit 
volume to obtain the density of cosmic ray sources. 
Figure 1 shows that about 1/3 of supernova explo- 
sion kinetic energy E goes to cosmic rays, which is 
in agreement with the empirical model of cosmic 
ray origin. 

The spectrum of accelerated energetic ions 
other than protons has the same shape if expressed 
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Fig. 1. — Source spectra produced by supernovae 
Type la (solid line), Type IIP (dash line), Type 
Ib/c (dotted line), and Type lib (dash-dot line) 
assuming that the accelerated particles are pro- 
tons. 

as a function of magnetic rigidity Q{p/Z) with the 
appropriate absolute normalization determined by 
the injection process at thermal energies. 

The relativistic ions released into interstellar 
space from numerous supernova remnants, diffuse 
in galactic magnetic fields, interact with interstel- 
lar gas, and finally escape through the cosmic ray 
halo boundaries into intergalactic space, where the 
density of cosmic rays is negligible. The main 
characteristics of cosmic ray propagation in the 
Galaxy needed for calculation of the cosmic ray 
spectrum in the so called leaky-box approxima- 
tion is the escape length X e , that is the average 
matter thickness traversed by cosmi c rays before 
they e xit from the Galaxy, see e.g. IStrong et al. 
([20071 ) . The cosmic ray intensity obeys the rela- 
tion / oc v sn Q{X~ l + a/m a ) , where a is the 
nuclear spallation cross section for a given type of 
relativistic nuclei moving through the interstellar 
gas, m a is the mean interstellar atom mass. 

The value of the escape length is determined 
from the relative abundance of secondary nuclei 
(primarily from the Boron-to- Carbon ratio ) incos- 
mic rays. Based on the paper Jones et al.l ( 2001 ). 
we choose the escape length in the form 



X e = 11.8 (v/c) (p/4.9ZGV) 



-0.54 



l/cm 2 (7) 



at p/Z > 4.9 GV; and X e oc v/c at p/Z < 4.9 
GV. Eq. (7) means that the resulting spectrum is 
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steeper than the source spectrum by 0.54 at high 
enough energies. We use Eq. (7) at all energies 
despite the fact that the B/C ratio is not measured 
above 100 - 1000 GeV. 

The source normalization for nuclei from pro- 
tons to Iron was made in our calculations from 
the fit to observed cosmic ray composition at one 
reference energy 10 3 GeV. 

The results of the calculations are shown in Fig- 
ure 2a for the interstellar spectra of protons, He- 
lium and Iron nuclei at kinetic energies per nucleon 
1 GeV/n< E < 10 3 GeV/n where the charge res- 
olution of cosmic ray experiments is high and the 
escape length X e (E) is known. The agreement 
between our theoretical predictions and the ob- 
servations of cosmic-ray energy spectra supports 
the validity of our acceleration and propagation 
models. (The discordance between the calculated 
interstellar and the measured at the Earth spec- 
tra at energies below 10 GeV/n is due to the solar 
wind modulation effect.) 

The combined spectrum of all protons and ions 
with particle energies E > 10 3 GeV is shown 
by the solid line in Figure 2b. It was assumed 
that the charge composition of accelerated parti- 
cles was the same in all types of SNRs except that 
the highest-energy part of the spectrum produced 
by Type Ib/c supernovae at pc/Z > 10 GeV had 
no hydrogen. This reflects the composition of pre- 
supernova, the Wolf-Rayet star with a fast H-poor 
wind. It was assumed that the shock produced 
by a supernova explosion first propagates through 
the 5 solar masses of the Wolf-Rayet star mate- 
rial uniformly spread in the central part of the 
wind bubble and propagates through the mate- 
rial with normal composition after that. It may 
explain why the protons do not do minate in the 



cosm i c-ray composition a t the knee ([Antoni et al 



20051 iBudnev et alJliooi ). In fact the proton-to- 
Helium ratio in cosmic rays is probably decreasing 
with energy st arting at s maller energies 10 2 — 10 3 
GeV/nucleon (|Seoll2007l) . We plan detailed inves- 
tigation of this point in a separate paper. 

The calculated spectra show remarkably good 
overall fit to observations up to energies ~ 5 • 10 9 
GeV, where the transition to extragalactic cos- 
mic rays with the characteristic GZK suppression 
above 3 x 10 10 GeV likely occurs. To a good ap- 
proximation the bending of the observed spectrum 
at around the knee energy 3 • 10 6 GeV is repro- 




E GeV/nucleon 




E, GeV/particle 

Fig. 2. — (a) The calculated interstellar (not 
corrected for solar modulation at low ener- 
gies) spectra of protons, Helium, and Iron 
below energy 10 3 GeV/n. Observational 
da ta from AM S 
2 (jPanov et al 



( Alcaraz et ail |200C|). ATIC 



2006 ), BESS-TeV ( Haino et al 



20041 ) . HEAO-3 (En gelmann et al.l 119901 ). and 
TRACER (|Ave et al.l 120081) experiments are 
shown, (b) The all particle spectrum above 10 3 
GeV calculated under the assumptions that escape 
length is determined by Eq. (7) at all energies 
(solid line) and has a cut o ff at 2 • 10 7 Z GeV (dash 
line) . See iBlumer et al. (120091) for references to 
the observational data shown by grey symbols. 
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Fig. 3. — Calculated mean logarithmic mass of 
produced in the Galaxy cosmic rays (thick solid 
line) compared to observational data based on the 
average depth of sho wer maximum as presented in 
iBlumer et all (|2009h . 



duced although no special efforts were made to 
force the theory to fit the data. The bending is 
due to the combined effect of the summation over 
different types of SNRs and over different types of 
accelerated nuclei. 

The complicated chemical composition of high 
energy cosmic rays is illustrated in Figure 3 where 
the calculated mean logarithmic atomic number 
of cosmic rays < \n(A) > is presented. The in- 
crease of < ln(.A) > at energies from fO 5 GeV to 
1 7 GeV is due to the dependence of the knee po- 
sition on charge Pknoc oc Z for each kind of ion 
accelerated in Types la, IIP, Ib/c SNRs. Type 
lib SNRs with normal composition dominate at 
rigidities p/Z > 5 • 10 6 GV. They have a knee at 
about Pk nee /Z rj 5 • 10 7 GV and provide progres- 
sively heavier composition to the very high ener- 
gies. It should be pointed out that the increase of 
< ln(A) > at E > 3 x I0 8 GeV predicted in our 
calculations is not supported by the available ob- 
servations. If confirmed, these observations may 
signify the dominant contribution of extragalactic 
cosmic rays with light composition at these ener- 
gies. 

The obtained cosmic ray spectrum shown by 
the solid line in Figure 2b is very attractive for 
the explanation of cosmic ray data. However, the 
use of the escape length (7) at ultra high energies 
is not justified. Experimentally, the value of X e is 
determined from the abundance of secondary nu- 
clei in cosmic rays with good statistics only up to 



about 100 GeV/n (jStrong et all 120071 ). If cosmic 
ray transport in the Galaxy is described as dif- 
fusion, the diffusion coefficient can be expressed 
through the escape length as D ss v[iH/2X e (here 
|U s» 0.003 g/cm 2 is the surface mass density of 
Galactic gas disk, H « 4 kpc is the height of 
the Galactic cosmic-ray halo), which gives D ps 
1.3 • 10 28 (pc/ZGeV) - 54 cm 2 /s. The diffusion ap- 
proximation can be used when the diffusion mean 
free path 3D/v is less than the size of the system 
if, which results in the condition pc/Z < 2 ■ 10 7 
GeV. At somewhat higher energies the particles 
accelerated in the galactic disk fly straight away 
from the Galaxy with a flat (source) energy spec- 
trum and close to hundred percent anisotropy. 
Certainly this picture does not represent the re- 
ality that may be due to the strong intermittency 
of very high energy cosmic rays produced by ran- 
dom short bursts of not very numerous sources. 
The dash line in Figure 2b shows the results of 
calculations made under the assumption that cos- 
mic ray particles with energies pc > 2 • 10 7 Z GeV 
freely escape from the Galaxy without being de- 
tected by observer at the Earth. The predicted 
spectrum may fit observations only below about 
5 • 10 7 GeV in this case. 

The validity of diffusion approximation extends 
to higher energies in the diffusion model with dis- 
tributed reacceleration on the interstellar turbu- 
lence where X* oc (p/ Z) -1 / 3 at high rigidities, see 
Seo fc Ptuskinl (|1994l ). However, this scaling does 



not reproduce the observed cosmic ray spectrum 
for the calculated source spectrum, see also dis- 
cussion in Section 3. It is worth noting that the 
uncertainty in our knowledge of parameters of the 
interstellar turbulence does not allow to decide be- 
tween two basic models of cosmic ray propagation: 
the plain diffusion model and the d iffusion model 
with di stributed reacceleration 
(|2006bh . 



see 



Ptuskin et al 



The physical pattern of cosmic ray propagation 
is different in the models with a Galactic wind. 
The wind model with selfconsistently calculated 
cosmic-ray transport coeffici ents reproduces we ll 
the data on secondary nuclei ( Ptuskin et al.fl997j ) . 
The supersonic wind is probably terminated by 
the shock at ~ 0.5 Mpc from the Galactic disk. 
The confinement of very high energy cosmic rays 
in the Galaxy can be more efficient in this model as 
compared to the diffusion model with a static flat 
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halo of the size ~ 4 kpc discussed above. In any 
case, trajectory calculations in galactic magnetic 
fields are needed to study cosmic ray propagation 
when the diffusion approximation breaks up at ul- 
tra high energies. The detailed consideration of 
this issue is beyond the scopes of the present pa- 
per. 

3. Discussion and Conclusion 

We have calculated the steady state spectrum 
of cosmic rays produced b y SNRs in the Galaxy. 
Our n ew numerical code ( Zirakashvili fc Ptuskinl 
l2009al fbh for modelling of particle acceleration by 
spherical shocks with the back reaction of cosmic 
ray pressure on the shock structure was used in 
the calculations. The significant magnetic field 
amplification in young SNRs inferred from the ob- 
servations of the ir synchrotron X-ray radiations 
(jVolk et al.l 120051 ). and most probably produced 
by cosmic ray streaming instability, was also intro- 
duced in the calculations. It led to the inclusion 
of the Alfvenic drift in the equation for particle 
transport downstream of the shock. Four different 
types of SNR s with relative burst rates taken from 
Smartt et all (|2009h were included in the calcula- 



tions. The escape length Eq.(7) from the work of 



Jones et al.l (|2001l) was used to describe the prop- 



agation of cosmic rays in the Galaxy in the leaky- 
box approximation. The normalization to the ob- 
served intensity and chemical composition of cos- 
mic rays was made at 10 3 GeV energy. 

The results are illustrated by the solid lines in 
Figures 2 when Eq. (7) for X e is used without lim- 
itation and by the dash line when it is limited by 
the applicability of the diffusion approximation for 
cosmic ray propagation in the diffusion model with 
a flat static halo. The solid lines reproduce well 
the entire cosmic ray spectrum up to ~ 3 • 10 9 GeV 
while the dash line makes it up to less than ~ 5 TO 7 
GeV. Further investigations of cosmic ray propa- 
gation in galactic magnetic fields at ultra high en- 
ergies are needed to refine the predicted shape of 
the spectrum produced by the Galactic SNRs at 
ultra-high energies. This is important in light of 
the discussion about transition from the Galactic 
to extragalactic co mponent in the observ e d cos- 
mic r ay spectrum (jBerezinskv et al. l2006t iHillasI 
20061) . 

Our results can be compared to the earlier work 



Berezhko fc Volkl (|2007l ). where the Alfvenic drift 
effect was not taken into account and only Type 
la SNRs that represents about 30% of all SNRs 
and can not efficiently accelerate particles to ultra- 
high energies were considered. The absence of 
Alfvenic drift led to a very flat cosmic ray source 
spectrum, which required too strong dependence 
of the escape length on rigidity X e oc (p/Z)~°' 75 
inconsistent with the data on secondary nuclei. 
In fact, this problem with a too hard predicted 
source spectrum was the main motivation for the 
present work. We included the effect of Alfvenic 
drift downstream of the shock and found that the 
resulting source spectrum fits the observations. To 
show the importance of this effect clearly, we ac- 
cepted basically the same set of other parameters 
and assumptions of the model of shock a ccelera- 
tion in SNRs as berezhko & Volkl (|2007h . Also, 
three additional types of SNRs were included into 
the consideration. 

The present work demonstrated that super- 
novae can in principle produce the source spec- 
trum of galactic cosmic rays required by the em- 
pirical model of cosmic ray origin. More work is 
needed to understand how robust are our results. 

First of all it is more detailed analysis of MHD 
effects that accompany the process of the dif- 
fusive shock acceleration in evolving SNRs. In 
particular, a further investigation of generation 
and transport of MHD turbulence at astrophys- 
ical shocks is necessary. For simplicity we de- 
scribed the properties of MHD turbulence in our 
calculations by only one parameter Ma- Its rep- 
resentative value Ma = 23 corresponding to the 
typical mag nitude of amplified magnetic field de- 
termined by Volk et al] (|2005f ) from the observa- 
tions of non-thermal X-ray emission from young 
SNRs gives the particle spectrum that is in accor- 
dance with cosmic-ray observations. However we 
checked that this result strongly depends on the 
assumed value of Ma which regulates the spectral 
slope via Alfvenic drift downstream of the shock. 
As an example, for Ma = 15 the spectrum of 
accelerated particles is so steep that the escape 
length X e cx (p/Z) -1 / 3 typical for the propaga- 
tion model with distributed reacceleration would 
reproduce the observed cosmic ray spectrum. The 
problem is however that the streaming instability 
of particles with steep spectrum can hardly gener- 
ate the strong magnetic field at large scales that 
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is needded for an acceleration of high-energy par- 
ticles. All this sets one thinking about possible 
self-consistent mechanism that maintains the spe- 
cific values of Ma- 

The Alfvenic drift effects downstream of the 
shock will be verified in the nearest future via com- 
bined observations of young supernova remnants 
in Fermi, HESS, MAGIC, VERITAS and some 
other experiments. The observations of older rem- 
nants will also help to understand how the max- 
imum energy of accelerated particles depends on 
the age of a remnant. Although the maximum 
energy pro bably quickly decreases in the old rem- 
nants (see Ptuskin fc Zirakashv ili (2003)), we ne- 
glect this effect in our present calculations since 
it would require an introduction of an additional 
rather uncertain parameter. 

A further work is also necessary for determi- 
nation of the dependence of injection efficiency ry 
on the shock parameters. The injection efficiency 
has strong influence on the resulting particle spec- 
tra. It is because the presence of the Alfvenic drift 
downstream of the shock steepens the spectrum 
of accelerated particles. So, it is more difficult 
to modify the shock by the pressure of high en- 
ergy particles compared to the case without the 
Alfvenic drift. 

The assumed strong heating in the shock pre- 
cursor also influences the resulting particle spec- 
trum. Although the main part of the cosmic ray 
energy that goes to the Alfven wave generation 
is eventually transformed into the gas heating, 
the level of Alfvenic turbulence {8B/B ~ 1) is 
high enough to provide the scattering of cosmic 
ray particles in the regime close to the Bohm dif- 
fusion. Note that quasi-linear estimates of the 
Alfven wave amplitu de give 6B/B >> 1 for yo ung 
supernova remnants ( McKenzie fc Volk |[l982 ). In 
reality Alfven waves generated by the resonant 
streaming instability probably stops when SB/B 
reaches the value close to 1 because of collisions 
of the moving adjacent plasma eleme nts similar 



the same type of supernovae are probably too sim- 
plified to correspond to reality. More compre- 
hensive analysis taking into account earlier work 



to results of 



The gas 



heating is then provided by weak shocks produced 
in these collisions. 

Our assumption that the composition of accel- 
erated particles is the same for all types of nu- 
clei (except the very high energy part of the spec- 
trum produced by Type Ib/c SNRs) and our ig- 
noring the dispersion of SNR parameters within 



nensivc analysis taxin g into account earlier work 
(Silb erberg et al.lll99ltlSveshnikoval2003l;|Popescu 



20071 ) is needed, although it requires the introduc- 
tion of additional not well-known astrophysical pa- 
rameters. 

Another difficult problem is accounting for the 
fluctuations resulting from the discrete nature 
of supernovae in space and time. Fluctuations 
of cosmic ray intensity a nd anisotropy were ad- 
dresse d in particular bv lErlvkin fc Wolfendalel 
(|2006l ): IPtuskin et al.l (|2006al ). 
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